arXiv:1502.01706v2 [cond-mat.quant-gas] 29 Apr 2015 


Strongly Interacting Quantum Gases in One-Dimensional Traps 


Li Yang^, Liming Guan^’^, and Han Pu^’^ 

^Department of Physics and Astronomy, and Rice Quantum Institute, Rice University, Houston, TX 77251, USA 
^Institute for Advanced Study, Tsinghua University, Beijing, 100084, P- R- China 
^ Center for Cold Atom Physics, Chinese Academy of Sciences, Wuhan 430071, P. R. China 

Under the second-order degenerate perturbation theory, we show that the physics of N particles 
with arbitrary spin confined in a one dimensional trap in the strongly interacting regime can be 
described by super-exchange interaction. An effective spin-chain Hamiltonian (non-translational- 
invariant Sutherland model) can be constructed from this procedure. For spin-1/2 particles, this 
model reduces to the non-translational-invariant Heisenberg model, where a transition between 
Heisenberg anti-ferromagnetic (AFM) and ferromagnetic (FM) states is expected to occur when the 
interaction strength is tuned from the strongly repulsive to the strongly attractive limit. We show 
that the FM and the AFM states can be distinguished in two different methods: the hrst is based 
on their distinct response to a spin-dependent magnetic gradient, and the second is based on their 
distinct momentum distribution. We confirm the validity of the spin-chain model by comparison 
with results obtained from several unbiased techniques. 
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I. INTRODUCTION 

One dimensional (ID) quantum systems have received 
much attention during the past many decades. This is 
due to the fact that quantum effects are more pronounced 
in reduced dimensions, and also to the fact that many 
ID models, such as the Lieb-Liniger model [l| and the 
Gaudin-Yang model SJ3, can be solved exactly with 
Bethe ansatz method [JjQ- Most exactly solvable models 
require the underlying systems to be translation invari¬ 
ant and the models can then be integrable. The presence 
of an external trapping potential in general breaks the 
integrability. One notable exception to this is a system 
of ID spinless bosons with infinite contact repulsion (the 
so called Tonks-Girardeau gas) confined in an arbitrary 
trapping potential, which can be mapped into a non¬ 
interacting spinless Fermi gas @,0 , and has been realized 
in experiments using ultracold atoms However, if 

the particles possess spin degrees of freedom, the prob¬ 
lem becomes much more complicated. In recent years, 
there has been works on constructing the ground state 
of ID spinful bosons and fermions with infinite or nearly 
infinite contact interaction [IMl- It has been shown 
that, at exactly infinite interaction, the ground state of 
such spinful particles possesses degeneracy as the energy 
is independent of the spin configuration. Slightly away 
from this infinite repulsion limit, a perturbation theory 
can be constructed using 1 /g (where g is the strength of 
the contact interaction) as the small parameter. In this 
way, the ground state is governed by an effective Hamil¬ 
tonian defined within this degenerate subspace Sin. 

In this work, we will explicitly construct an effective 
model for ID strongly interacting particles using a per¬ 
turbation approach. Here the unperturbed system con¬ 
sists of particles with infinite contact interaction, i.e., 
g = oo. At finite but large |( 7 |, we take 1/g as the 
small perturbation parameter. We will show that we 
need to take the perturbation to second order in order to 


break all the spin degeneracy. In this way, we can con¬ 
struct an effective Hamiltonian which takes a form of a 
non-translational-invariant Sutherland model [l8| , which 
arises from the effective super-exchange interaction be¬ 
tween neighboring particles. One can intuitively under¬ 
stand the emerence of the super-exchange term as follows. 
At g = oo, particles are inpenetrable in ID and they can¬ 
not exchange positions with their neighbors. Away from 
g = oo, there will be small but finite probabiliy that 
two neighboring particles can exchange positions, which 
gives rise to the effective super-exchange interaction. For 
spin-1/2 fermions, which we focus on in this work, the 
exchange operator can be written in terms of spin oper¬ 
ators, and the Sutherland model reduces to the Heisen¬ 
berg model. It immediately follows that the ground state 
of spin-1/2 fermions is a Heisenberg anti-ferromagnetic 
(AFM) state in the strongly repulsive limit, and a fer¬ 
romagnetic (FM) state in the strongly attractive limit 
(we exclude the tightly bound molecular states on the 
attractive side, i.e., we consider the upper branch of the 
system). We investigate the properties of such a system 
and demonstrate experimental signatures that allow us to 
distinguish the AFM and the FM states, using both the 
effective model and several unbiased methods, and show 
that the former is indeed valid in the strongly interacting 
regime. 

The main advantages of the effective model are two 
fold. First, from a conceptual point of view, the effective 
model provides new insights to the quantum magnetic 
properties of strongly interacting particles in ID. Sec¬ 
ond, from a practical point of view, the effective model is 
much easier to handle in comparison to unbiased meth¬ 
ods. As a result, the effective model allows us to deal with 
more particle numbers and to investigate the dynamics 
to longer time scales. To this end, we benchmark our ef¬ 
fective model against several unbiased methods and show 
that the former is indeed valid in the strongly interacting 
regime. These benchmark calculations also demonstrate 
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that calculations based on the effective model are much 
more efficient and take much less time than those based 
on unbiased methods. 

The rest of the paper is organized as follows. In 
Sec. II, we derive the effective spin-chain Hamiltonian us¬ 
ing a second order perturbation theory. We compare the 
energy spectrum obtained from this Hamiltonian with 
that obtained from a numerically exact Green’s function 
method. In Sec. HI, we calculate the density profiles 
of the ID trapped system in both real and momentum 
spaces. We show that the FM and the AFM states pos¬ 
sess identical real space density profile, but with distinc¬ 
tive momentum distribution. In Sec. IV, we study the 
system’s response to a spin-dependent magnetic gradient, 
which breaks the SU(2) symmetry and hence mixes the 
AFM and the FM states. In Sec. V, we show how the spin 
symmetry breaking term helps to realize the FM state in 
practice. Finally, in Sec. VI, we discuss the advantages of 
the effective model over those unbiased methods, which 
serve as an important motivation for this work. Many of 
the technical details can be found in the Appendices. 


II. EFFECTIVE SPIN-CHAIN MODEL 


We consider a one-dimensional system with N strongly 
interacting spinful particles with mass m trapped in an 
arbitrary external potential, with the Hamiltonian 


N 

i=l 

Hf 


1 
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Here we have set h = m = 1. For infinite repulsion 
the particles become impenetrable and behave like spin¬ 
less fermions. If the N particles are spinless bosons, the 
many-body wave function can be constructed by Bose- 
Fermi mapping [^. For spinful fermions, the correspond¬ 
ing wave function can be generalized |ll| to 


■ • •XAr,ai ■ ■ -fTAr) = [^PaO^ 0 x] , (2) 

P 


where (pA is a Slater determinant which represents the 
eigen-wave function of N spinless fermions governed by 
Hamniltonian Hf. Here 9^ is a sector function (i.e., gen¬ 
eralized Heaviside step function) of spatial coordinates, 
whose value is one in spatial sector Xi < X 2 < ■ ■ ■ < Xn, 
and zero in any other spatial sectors, y is a spin wave 
function, and P is the permutation operator whose con¬ 
vention of acting on spatial and spin wave functions is 
presented in Appendix [Al 

To obtain an effective Hamiltonian for spinful fermions 
in the strongly interacting regime, we use the perturba¬ 
tion theory. To this end, we consider Hf as the pertur¬ 
bation, and iLint as unperturbed Hamiltonian. This is 
in the same spirit as the procedure for constructing the 
effective spin model from the Hubbard model in the large 


interaction limit [l9j | . The unperturbed Hamiltonian Hint 
has a degenerate ground state subspace with zero eigen- 
energy = 0. This subspace is the space of all the 
anti-symmetric wave functions satisfying the boundary 
condition 'i’xi=xj = 0 @[ni- Equation ([2]) with a full 
set of ipA^s constitute a complete basis for this subspace. 
We define a projection operator Vq into this subspace 
and its complementary operator Vi = 1 — Vq. Now let 
us consider the effect of Hf on this subspace under the 
framework of degenerate perturbation theory. The first 
order effective Hamiltonian reads H^^^ = VoHfVo- The 
ground states of H^^^ still form a degenerate subspace 
whose eigen-vectors take the same form as Eq. © with 
(fiA representing the lowest-energy Slater determinant for 
Hf. (Erom now on, we denote ipA as such a lowest-energy 
Slater determinant.) To lift the remaining spin degener¬ 
acy, we therefore have to carry out the perturbation cal¬ 
culation to second order. Let Qq be the projection op¬ 
erator into the ground state subspace of H^^\ Applying 
standard degenerate perturbation theory, we obtain the 
second-order effective Hamiltonian as (see Appendix [B] 
for details). 

H-G) ^ - ViHfQo. (3) 

-^int ~ ^int 

After some algebra (for details, see Appendix[C|), we find 
that, after neglecting a constant QoHfQo, the effective 
second-order Hamiltonian can be written as 

N-l 

HeS = — Ci{l — £i i+i), (4) 

9 U 

where £if+i is the exchange operator acting on a spin 
state X within the subspace defined by Qq, and its effect 
is to exchange the A and {i -I- 1)*'' particles, and the 
coefficients 

Ci = N\ f Y[ dxj \diipA\^ 5{xi+i - Xi)e\^^t,i ]: (5) 

3 

are positive constants independent of spin, where 
=9^l6{xi-Xi-i), 

is a reduced sector function (see Appendix ICl) . iLeff takes 
the form of the non-translational-invariant Sutherland 
model, and physically arises from the effective super¬ 
exchange interaction when g deviates away from infin¬ 
ity, as we have mentioned earlier. In the case of spinful 
bosons, following the same procedure leads an effective 
Hamiltonian similar to (|4|) with the minus sign before 
£ip+i replaced by the plus sign. This spin-chain model 
preserves the SU(2s-|-l) symmetry, where a single parti¬ 
cle has spin s. Being a bipartite Hamiltonian, the Lieb- 
Mattis theory is also satisfied. Since it is made up 

of permutation operators, it can also be block diagonal¬ 
ized in the irreducible representation of the permutation 
group Sn ho, mi. 







3 


We comment here that Eq. @ can be written in a 
different form, 'I' = with C{cr},p = 

l/(iV^! 7 V 4 ,!) ({cr}|-Px) being weights in different sectors 
for a spin configuration {a}. These weights can be 
regarded as variational parameters and determined by 
dE/dcp = 0 together with the Bethe-Peierls boundary 
condition [l^, . For strong but finite interaction, the 

eigen-energies read E = Eq — K/g + 0{l/g‘^), here K 
is the Tan contact [2^. An effective spin model can be 
constructed from this variational approach, as done by 
several groups [H, [13, [H, . Our result based on the 

perturbation calculation is consistent with these results. 

To benchmark the spin-chain model, we show in 
Fig.[TJa) the low energy spectrum of a three-body system. 
Similar benchmarks were also performed in Refs. |lu.[^. 
In this work, we focus on spin-1/2 fermions, and label 
the two spin species as f and The external potential 
is chosen to be a harmonic potential with frequency w. 
In our calculation, we take a; = 1 along with h and m, 
and the observables are normalized to dimensionless val¬ 
ues: X ~ x/yjhj{mu}), p ~ p/y/hmuj, and E ^ E/hu). 
The main figure of Fig. [TJa) is obtained by the unbiased 
Green’s function method based on the original many- 
body Hamiltonian o M- In the inset, we compare this 
exact spectrum (dots) with the spectrum obtained from 
the spin-chain Hamiltonian Hgg (solid lines). As one can 
see, in the strong interaction regime with \/\g\ 1 , the 

spin-chain model faithfully reproduces the exact spec¬ 
trum of the upper branch when the tightly bound molec¬ 
ular states on the attractive {g < 0 ) side are ignored. 

We can gain some insights into the spectrum of i?eff 
by noting that the eigenvalues of the exchange operator 
£i^i+i are ±1. Therefore, for 5 > 0, the spectrum of ideff 
has a lower bound of — ( 2 /g) (corresponding to 

a fully anti-symmetric spin configuration with £i^i+\ = 
— I for any i), and an upper bound of 0 (corresponding 
to a fully symmetric spin configuration with £i^i+i = I 
for any i). We remark that the fully anti-symmetric spin 
configuration can only be realized for 2s-|-1 > N. For not 
too small N, this requires a fermionic species with large 
spin s. Recent cold atom experiments have witnessed 
realization of high spin Fermi gases in alkali-earth atoms 
(26l - l^ . For g < 0, the spectrum is inverted and bound 
between 0 and \2/g\ Ci- 


III. DENSITY PROFILES IN REAL AND 
MOMENTUM SPACES 



l/f7 



FIG. 1: (Color online) Energy spectrum of the relative motion 
as a function of 1/g for three fermions with (A]-, N{) = (1, 2), 
without (a) and with (b) the spin-dependent magnetic gradi¬ 
ent. For (b), we have G = 0.05. The main figures are obtained 
using the Green’s function method. The red dotted lines in 
negative g area represent the tightly bound molecular states. 
The inset figures show the comparison between the spectrum 
obtained from the Green’s function method (dots) and that 
from the effective spin-chain model (solid lines) near 1/g = Q. 
In all the figures presented in this paper, we have adopted the 
trap units with ft = m = = 1. Consequently, the energy E 

is in units of hw, and the interaction strength g is in units of 
y/Wuj/m. 


where Bi are the Pauli spin matrices for the fth atom. 
Hence we can rewrite the effective Hamiltonian (l4|) as 

1 ^-1 

HeS = — Ci{l - Bi ■ dj+i)/2, (6) 

3 U 


Let us now examine in detail the density profiles in 
both real and momentum spaces for the ground state of 
ideff- For spin-1/2 fermions, the exchange operator can 
be written in terms of the spin operators: 

£i^j — (1 -j- CTj ■ (7j^l‘2 , 


which takes the form of the non-translational-invariant 
Heisenberg model with Ci/{2g) plays the role of the su¬ 
per exchange coefficient between the zth and the (z-l-l)th 
spin. The effective spin-spin interaction is ferromagnetic 
for g < 0, and anti-ferromagnetic for g > Q. We there¬ 
fore label the corresponding ground state FM for g < 0 
and AFM for g > 0, as shown in the inset of Fig. [TJa), 
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which is consistent with the Bethe ansatz result for the spin lowering operator, 
homogeneous case [13, HH . Note that as the number of 

atoms in each spin species are individually conserved, the To find the density profiles in both real and momentum 
spin configuration for the FM state here can be written spaces, let us first introduce the one-body density matrix 
as Itt • • ■ t), with S~ = being the total element defined as 

I 


Pa'a[x',x)= ^ j dX2-■ ■ dXN ^*{x\x2-■ ■ xn,cf', 02- ■ ■ (rN)^{x,X2 ■ ■ ■Xn,CF,(J2-■ ■ CFn) , (7) 

C72---ITN 


from which the real-space and momentum space density 
profiles can be calculated as 

Pa{x) = Npc^a{x,x), 

Pa{p) = {N/2 tt) J dx J dx'Pa,a{x\x). 

In Appendix [Dl we provide the details of calculating 
the one-body density matrix element given a many-body 
wave function as in Eq. 

In Fig. [21 we present the density profiles for = 2 
spin-1/2 fermions with (A^-|-,A^ 4 ,) = (1,1). For this two- 
body problem, exact analytic solutions for arbitrary in¬ 
teracton strength g can be found ■ Results desplayed 
in Fig. 121 are obtained from the exact method. As a re¬ 
sult, we are not limited to large |g|. Note that the FM 
state corresponds to a fully symmetric spin configuration 
X, and its density profiles, which are ^-independent, are 
identical to a system of N spinless fermions. More specifi¬ 
cally, Pa{x) = (Na/N) where (j)i{x) is the 

ith eigen-wave function of the single particle Hamilto¬ 
nian; and Pcr(p) decays as exp(—in the large p limit. 

The AFM state, on the other hand, possesses a fully 
anti-symmetric spin configuration and its density pro¬ 
files are sensitive to the value of p. As l/p —>■ 0, the 
real-space density profile of the AFM state approaches 
that of the FM state, whereas the momentum space den¬ 
sity profile remains distinct for these two states. Hence, 
in the strongly interaction limit, the density profiles for 
the AFM and the FM states are indistinguishable in real 
space, but distinguishable in momentum space. This 
statement remains true for N > 2. 

As a further example, we consider a system of 
(iV'|-,iVi,) = (4,4) spin-1/2 fermions in the strongly in¬ 
teracting limit. In Fig. Owe show the momentum space 
density profiles. The black dashed line corresponds to 
the momentum distribution of the FM state (which is 
the same as the momentum distribution of N spinless 
fermions), and the red solid line to that of the AFM state. 
The AFM state has a nonzero Tan contact K, and in the 
large momentum limit, we have p{p) = K/{2 ttp‘^) [33| . 
This is confirmed by our numerics as shown in the inset 
of Fig. 0 For comparison, we also show the momentum 
distribution of a fully anti-symmetric spin state, which 


AFM fm 

1 /g =-1 1/g=-0.1 l/g=0 1/g=0.1 1/g=1 



FIG. 2: (Color online) Real space density profiles (upper 
panel) and momentum space density profiles (lower panel) 
for (Af, Ni) = (1, 1) versus 1/g for AFM (red solid lines) and 
FM (black dotted lines) states. In our trap units, position x 
is units of \/hf{nujj), the real space density p{x) is in units 
of y/rnuTJh, the momentum p is in units of \/hmuj, and the 
momentum space density p{p) is in units of 


coincides with the momentum distribution of N spinless 
bosons in the Tonks-Girardeau limit. As we mentioned 
earlier, the fully anti-symmetric spin state is only possi¬ 
ble when 2s +1 > N [^. We emphasize again that these 
different states have identical real space density profile, 
but can be distinguished from their distinctive momen¬ 
tum distribution. 


IV. RESPONSE TO SPIN-DEPENDENT 
MAGNETIC GRADIENT 

The form of the spin-chain effective Hamiltonian i7eff 
makes it clear that a quantum phase transition is induced 
as l/g is tuned across zero, which can be achieved using 
the technique of confinement induced resonance (ssj, l33| . 
In practice, however, more effort is required to observe 
this phase transition. The AFM ground state for 5 > 0 
can be straightforwardly prepared. Such is not the case 
for the FM state on the attractive side with g < 0. This 
is due to the fact that, for g < 0 , there exist many bound 
molecular states with lower energies than the FM state, 
as can be seen from Fig. [TJa). If one simply prepare 
the system on the attractive side, these molecular states. 
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FIG. 3: (Color online) Momentum distribution for (A^-|-, Ni) = 
(4,4). The black dashed curve is for the fully spin symmet¬ 
ric FM state, which has the same momentum distribution 
as N spinless fermions. The red solid curve is for the AFM 
state. The blue dash-dotted curve is for the fully spin anti¬ 
symmetric state, which has the same momentum distribu¬ 
tion as N spinless Tonks-Girardeau bosons. The inset shows 
the momentum distribution for the AFM state in the large 
momentum limit, in comparison to the theoretical prediction 
Kf{2'Kp‘^) (green dotted curve), where K is Tan contact. 


not the FM state, will be realized. Hence to create the 
FM state, one needs to start from the AFM state on 
the repulsive side and adiabatically tune the interaction 
strength to the attractive side. However, the spin states 
are protected by symmetry: If we start from the AFM 
state and tune 1/g across zero, the system will remain 
as an AFM state and realize a fermionic super-Tonks- 
Girardeau state [ 13 , there is no coupling between 

the AFM and the FM states. To overcome this problem, 
we need to add a spin symmetry breaking term. One 
possibility is to add a spin-dependent gradient term. We 
will consider in detail how to realize the FM state in the 
next section. Here we first investigate how the AFM and 
the FM states respond to such a gradient term. 

To this end, we introduce a weak spin-dependent mag¬ 
netic gradient which adds a term — to the 

Hamiltonian o, where G, which we will take to be non¬ 
negative, chracterizes the magnitude of the magnetic gra¬ 
dient. The effective spin-chain Hamiltonian will be mod¬ 
ified corresondingly as 

^ N-l N 

HeS = X] ■ '^i+i)/2 - G'^Dia- , (8) 

^ i=l i=l 

where Di = N\ J Xi\(pA\‘^&^ H^i dxj represents the po¬ 
sition of the fth atom. In Fig. [Ijb), we plot the en¬ 
ergy spectrum for a three particle system in the presence 
of weak spin gradient, obtained from both the Green’s 
function method and the effective model. Again we see 
excellent agreement in the strongly interacting regime. 


Comparing the insets of Fig.[T](a) and (b), one can easily 
see that the gradient term lifts the spin degeneracy at 
1/5 = 0 , and the ground state is now separated from ex¬ 
cited states by a finite gap, which facilitates the adiabatic 
preparation of the FM state to be discussed later. 

The spin gradient tends to separate the two spin 
species (39|. To quantify this effect, we define 

1 ^ 

A = ^ ^{xiat ), (9) 

i=l 


which measures the center-of-mass separation between 
the two spin species. Here the expectation value is taken 
with respect to the ground state of the effective Hamilto¬ 
nian dH). In the absence of the gradient (G = 0), A = 0 
for both the FM and the AFM states. Under the effective 
spin-chain model, A is a function of Gg only. 

As a first example, we again consider a two particle 
system with = (1,1). For this simple system, 

Hamiltonian ([ 8 ]) can be easily diagonalized, and A has 
an analytic expression: 


A = 


2G|5| + x/l+4(G|5|)2 


- 1 


2G|5| +v'l + 4(G|5|)2 


Note that since A only depends on we conclude that 
the FM and the AFM state respond identically to the 
gradient in the two-body case. We plot A as a function 
of GI 5 I in Fig.|4|a). In the figure, we also plot the result 
obtained from an exact solution using the Green’s func¬ 
tion method with g = ± 20 , which are in good agreement 
with the effective model. The details of this solution can 
be found in Appendix [E] 

By contrast, for N > 2, the ground states for 5 > 0 
and g < 0 will response differently to the gradient. In 
Fig. [41(b) and (c), we plot A as a function of Gj^j for 
the cases (A^,A^) = (1,2) and (2,2), respectively. The 
dashed and solid curves correspond to the ground state 
of negative and positive 5 , respectively. In general, the 
ground state on the attractive side will have a stronger 
response. To benchmark the effective model, we stud¬ 
ied this problem using the Time-Evolving Block Decima¬ 
tion (TEBD) method [4?)l - l4^ . In TEBD, a many-body 
wave function is represented by a Matrix-Product state 
(MPS), which approximates a many-body wave function 
by making a truncation of the entanglement spectrum. 
Eor ID gapped system, whose entanglement is short- 
ranged, the truncation error is well controlled, and the 
TEBD method therefore represents an unbiased method 
and has been implemented widely to study ID systems. 
The symbols in Eig. |4l(c) are the TEBD results for posi¬ 
tive g. One can see that for large g, the results obtained 
from TEBD and the effective model agree with each other 
very well. 

To further quanitify the response to the gradient and 
show the difference between the AFM and the FM 
states, we define the magnetic gradient susceptibility as 
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l^((iA/dG)G=o, and the following relation can be readily 
derived: 


1 dA 
Ji\dG 


G=0 


2 

1 ^ 




E- 

n^O 


En — E() 


( 10 ) 


where \n) represents the nth eigenstate of the spin-chain 
Hamiltonian with G = 0, and En is the corresponding 
eigen-energy. |0) represents the ground state, which is the 
AFM (FM) state for positive (negative) g. In Fig. lUd) 
we plot this susceptibility as a function of the total par¬ 
ticle number N for the case with iV^ = = N/2. One 

can see that, as long as N > 2, the FM state possesses a 
larger susceptibility, i.e., is more prone to spin segrega¬ 
tion under the gradient, than the AFM state. Further¬ 
more, the susceptibility for the FM state grows rather 
rapidly as N increases, whereas that for the AFM state 
is not very sensitive to N. 




FIG. 4: (Color online) Separation between the two spin 
species as a function of G\g\ for (a) {N-f,Ni) = (1,1), (b) 
{Nf,Ni) = (1,2), (c) (Af^f) = (2,2). The black dashed 
curves are for the ground state with negative g, and the red 
solid curves are for the ground state with positive g. The sym¬ 
bols in (a) are obtained from the analytic solution detailed in 
Appendix IeI The symbols in (c) are TEBD results, (d) The 

susceptibility | as functions of N for Nf = Aj, = N/2. 

In our trap units, A is in units of \/ft/(ma;), and Gg in units 


V. ADIABATIC PREPARATION OF 
FERROMAGNETIC STATE 

In the previous section, we suggested a method of 
applying weak spin-dependent magnetic gradient to ap¬ 
proach the FM state in experiment. Here we will discuss 
the method in detail. The experimental protocol is in 
the following: (1) The system is initially prepared in the 


ground state with strong repulsion {g > 0) and a rela¬ 
tively large magnetic gradient. In the example presented 
in Fig. [5l we choose the initial values 1/g = 0.01 and 
G = 0.1. (2) From t = 0 to Ti, G is fixed at the initial 
value while the interaction strength is tuned to the at¬ 
tractive side as l/g{t) = 0.01 cos(7rt/Ti), which can be 
achieved with confinement-induced-resonance method. 
(3) Finally, from t = Ti to Ti + T 2 , g is fixed at its value 
at Ti, while the gradient strength G is slowly turned off. 
We vary G such that the instantaneous spin separation 
A follows the form 


A{t) = A[G(t)] 


A(Ti) cos^ 


7r(t - Ti) 

2T2 


( 11 ) 


The experimentally controlled parameters are plotted in 
Fig. [SJa) for Ti = 20 Tho and T 2 = 280 Tho, where Tho = 
2'k jixi is the harmonic trap period. 

In Fig.jSJb) we display the evolution of the spin separa¬ 
tion parameter A in an example system with (A|, N = 
(2,2), Ti = 20Tho, and Ti + T 2 = 100T^o, 200T^o, 
and 300 Tho. The dashed curves represent the target- 
ted instantaneous value as shown in Eq. (EH); while the 
solid curves are obtained by solving the time-dependent 
Schrodinger equation under the effective Hamiltonian 
Heff- As expected, the larger the T 2 , the better agree¬ 
ment between the solid and dashed curves. In the inset, 
we also show the fidelity, which is the overlap between the 
calculated wave function from evolving the Schrodinger 
equation and the instantaneous ground state wave func¬ 
tion given the values of g and G at the moment, for the 
case Ti -I- T 2 = 300 Tho. One can see that an FM state 
can be realized with very high fidelity. For a shorter total 
evolution time with Ti -|- T 2 = 100 Tho, we still obtain a 
hdelity higher than 94%. 

Although we have proposed to use a spin-depedent 
magnetic gradient to break the spin symmetry and fa¬ 
cilitate the adiabatic preparation of the FM state, in re¬ 
ality any spin symmetry breaking term can do the job. 
Experimentally, this mean one needs to introduce some 
perturbation to the system to which the two atomic spin 
states will respond differently. A possibility is to apply 
an off-resonant light with proper polarization such that it 
induces different light shift to different atomic spin states. 
This idea has been recently implemented to create spin- 
dependent optical lattices for cold atoms [ 13 , . 

Finally, we comment on the stability of the FM state. 
Due to presence of the tightly bound molecular states 
on the attractive side, the FM can only be metastable. 
In 2009, Haller el al. realized such a metastable state 
in a system of spinless bosons and the resulting 

state is the so called super Tonks-Girardeau (sTG) gas. 
In that experiment, a typical lifetime of about 100 ms 
is found. We expect the lifetime of the FM state in a 
spin-1/2 Fermi gas should be longer than the bosonic 
sTG gas. This is because the low-lying molecular states 
for fermions must be spin singlet. Therefore the spin 
symmetric FM state will be protected by spin symmetry 
against decaying into the molecular states. 
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FIG. 5: (Color online) Adiabatic preparation of the FM state. 
At t = 0, the system is prepared in the ground state with 
Xjg = 0.01 and G = 0.1. (a) The value of experimentally 
controlled parameters l/g(t) and G(t) for a total adiabatic 
evolution time SOOTho. (b) The solid lines represent A(t) ob¬ 
tained by solving the time-dependent Schrodinger equation 
under the effective Hamiltonian iileff. The dashed lines repre¬ 
sent Eq. which is the A(t) of the instantaneous ground 
state for the given values of g{t) and G{t). Three different 
total adiabatic evolution time is calculated, 100ThoiSOOTho, 
and 300 Tho. The inset shows the hdelity of the adiabatically 
prepared state for the total evolution time 300 Tho. 


VI. DISCUSSION 

We have shown here, for large interaction strength 
the original Hamiltonian Eq. m can be map into a spin- 
chain model governed the by the effective Hamiltonian 
Hes in the form of Eq. which is expected to com¬ 
pletely describe the physics of the upper branch in the 
strongly interacting regime. The great advantage of the 
effective model is that (1) it provides valuable insights 
into the quantum magnetic properties of strongly inter¬ 
acting one dimensional quantum gases, and (2) it is much 
easier and more efficient to solve in comparison to the 
original many-body Hamiltonian. We have benchmarked 
the static properties of the effective model with several 
unbiased methods (see Fig. [T] and Fig.|T|). 

As we have mentioned earlier, recently several other 
groups have obtained the same spin-chain effective 
Hamiltonian using a variational method [Il,[i3[2i,[23. 
Our perturbational approach is inspired by the similar 
technique used to construct effective spin models from 
Hubbard Hamiltonian in the large-C/ limit. Using this 
technique, the super-exchange interaction arises natu¬ 
rally. The Hubbard Hamiltonian describes lattice sys¬ 
tems. Our work thus broadens this approach to a con¬ 
tinuum model. From the perturbation calculation pre¬ 
sented in this work, we may readily obtain many-body 


wave functions accurate to order 1/g. Furthermore, it is 
in principle possible to extend the perturbation approach 
to higher orders to obtain more accurate results. These 
features will be exploited in the future to study more 
detailed properties of the system. 

In Fig. [5] we present another example. Here we con¬ 
sider a quench dynamics in which the system is ini¬ 
tially prepared in the ground state with 1/g = 0.01 and 
G = 0.05. At t = 0, the spin gradient is suddenly turned 
off and the evolution of the center-of-mass separation be¬ 
tween the two spin species A is calculated by solving 
the time-dependent Schrodinger equation. We solve the 
Schrodinger equaton using both the effective spin-chain 
model governed by HeS, and the TEBD method governed 
by the original many-body Hamiltonian. As can be seen 
from Fig. [51 the effective model nicely reproduces the 
TEBD result. We therefore demonstrated that the spin- 
chain model can be applied to study the dynamics of the 
system. This example also serves to showcase the advan¬ 
tages of the effective model in the dynamical situation: 
due to its smaller Hilbert space, it can capture much 
longer time scale behavior of the system. Furthermore, 
it takes a few days to obtain the TEBD result as dis¬ 
played in Fig. [6l in comparison to a few tens of seconds 
for the spin-chain result. 



tig 

FIG. 6: (Color online) Quench dynamics for {N-f, N^) = (2, 2). 
The initial state is prepared as the ground state with 1/g = 
0.01 and G = 0.05. At t = 0, G is set to zero, and the system 
starts to evolve in time. The red solid line is A(t) calculated 
using the TEBD method, and the blue dashed line is cal¬ 
culated using the effective spin-chain model. The inset figure 
shows the evolution for much longer time under the spin-chain 
model. In our trap units, t/g is in units of ^Jm/{h^uJ^). 

Another great advantage of the spin-chain model is its 
wide applicability M- The effective Hamiltonian o is 
valid for spinful fermions, and by changing the minus 
sign in front of the exchange operator it describes 

strongly interacting bosons. The coefficients Ci, as given 
in Eq. only depend on the total number of atoms and 
the external trapping potential, and are independent of 
whether the particles are bosons or fermions, nor are they 
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FIG. 7: (Color online) Dimensionless coefficients Ci for N = 8 
and 13, calculated using the Monte Carlo integral method 
(Veges Algorithm [s^). The solid lines are obtained using 
the approximate expression US- 


dependent on the single particle spin s. The formalism to 
derive the effective Hamiltonian is independent of particle 
numbers N. Hence it works for any N. However, for N 
particles, each coefficient Ci invovles an iV-dimensional 
integral, which becomes quite difficult to evaulate as N 
increases. In Ref. [1^, the authors conjectured that, for 
a harmonic trap, these coefficients are given by 


N{N-l)/2 


( 12 ) 


where K is the Tan contact for the AFM state corre¬ 
sponding to N^) = (TV — 1,1). In Fig. [71 we plot the 
calculated Ci for TV = 8 and 13 (symbols), in comparison 
with the above expression (lines), and find good agree¬ 
ment. Hence, at least for harmonic trapped systems, 
once we know the Tan contact, all the Ci coefficients can 
be obtained approximately using Eq. (HU. We should 
also remark that recent experimental progress has made 
it possible to investigate few-particle cold atom systems 
with well controlled particle numbers in the lab [4^ [^ . 

Finally we comment that we have considered here a 
system of ID trapped spinful particles with strong con¬ 
tact interaction, and assumed that the interaction is spin- 
independent (i.e., SU(2s-|-l) symmetric), characterized 
by a single interaction parameter g. It is possible, within 
the framework of the perturbation method developed 
here, to generalize the formalism into a situation with 
spin-dependent interaction strengths, as long as all in¬ 
teraction strengths are sufficiently large [l^. Finally, it 
is also possible to generalize our work to Bose-Fermi mix¬ 
tures [3 [^ , which can be compared with recent few- 
body studies of such mixtures [SOi - l^ . We will consider 
these generalizations in a future work. 
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discussions. This work was supported by NSF and the 
Welch Foundation (Grant No. C-1669). L.G. acknowl¬ 
edges support from the Tsinghua University Initiative 
Scientific Research Program. 


Appendix A: Convention for the permutation 
operators 


This section contains the convention about the permu¬ 
tation operators and its action on spatial and spin wave 
functions. A permutation operator P can be expressed 
as 


f 1 2 ••• TV \ 

{ Pi P2 ■■■ Pn ) ’ 


(Al) 


which means that the original particle index T, after the 
permutation, is changed into Pi. 

The action of the permutation operator P on a spatial 
wave function is defined by 


P'iIj{xi,X2 - ■ ■ xn) = xp^). (A2) 


Similarly its action on a spin wave function is defined by 


Px{<xi,(X2---(Tn) = x(c^Pi,o-P 2 ■■■<xpn), (A3) 


where Ui stands for the spin state for Tth particle. The 
spin wave function x is a rank-TV SU(ri) tensor with n = 
2s -I-1, if all the particles are spin-s particles. 

A general spin state can be written as superpositon of 
basis tensors (or spin Fock states). A basis tensor can 
be written as = I/1/2 ■ • • /w), which 

means the Tth spin is in fi state. By definition, the per¬ 
mutation operator acting on a spin basis yields 

P 1/1/2 ■■■ In) = PSaifJa^f^ ■ ■ ■ 

~ fl ^<^P2 h ' ' ' ^‘PPn 
~ 1 Sa2f 1 • ■ ■ -1 ^ ^ 

= l/pyi/p-1 ■ 

We denote £ij as the exchange permutation operator, 
which simply exchanges indices i j'. 

£z ,3 1/1 ••• /i ••• /j ••• /at) = 1/1 ••• /j •••/*••■ /w) ■ 

(AS) 

We also denote the symbol {m - ■ -n) as a loop permu¬ 
tation operator, which, if m > n (to < n), permutes 
the indices by to p- to-|-1 p- m + 2 ■■ ■ •<— n P- m 
(to •<—to—!•<—to — 2----<—n-<— to). 


Appendix B: Second-order degenerate perturbation 
theory 

Consider a Hamiltonian 

H = Ho+ V, (HI) 
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where Hq is the unperturbed Hamiltonian, which possess 
a degenerate manifold with eigen-energy We define 
Vq the projection operator onto this degenerate subspace. 
V represents the perturbation Hamiltonian. To calculate 
the zeroth-order wave function and first-order energy cor¬ 
rection, we need to diagnalize VqVVo in Vo subspace. 
Suppose that the first-order energy spectrum still con¬ 
tains a degenerate manifold with energy and we 

define Qq as the projection operator onto this remain¬ 
ing degenerate subspace (obviously Qq C Vq)- To lift 
the degeneracy in Qq, we have to consider second-order 
perturbation. 

To distinguish the states in Qq by energy, we need to 
diagonalize the following operator in the Qq subspace. 


By doing this, we can obtain the zeroth-order wave func¬ 
tions and second-order energy correction E^^\ To 
calculate the first-order wave function correction 
we can use the following two formulas: 


(B4) 


Qom-^(5y^PiFQo, (B2) 


where Vi 


I — Vq is the complimentary space to Po [54 1. 


where Qi = Vq — Qq is the complimentary space of Qo 
in Vq. The first-order wave function correcton within the 
Qo subspace can be fixed to be zero, because we have a 
freedom before normalize the total wave function. 


I 

Appendix C: Derivation of the effective spin-chain model 


For N particles with contact interaction in a trap, the Hamiltonian is given in the main text as Eq. O: 


N r 




i=l 


1 52 
' 2 ^ 


+ V{xi) 


Hf 


g^5{xi-Xj) . 
i<j 

^ ^ > 

Hint 


(Cl) 


In the strongly interacting regime, we take the interaction Hamiltonian Hint as the unperturbed Hamiltonian, and 
the single-particle Hamiltonian Hf as perturbation. The ground state Hint is degenerate with energy E^^^ = 0. 

This section details how to derive the second-order perturbation effective Hamiltonian (IB2I) . with Hq replaced by 
Hint and V replaced hy Hf, into a spin-chain model. First consider the operator ViHfQo acting on an arbitrary state 
in Qo in the form of Eq. © in the main text, which we label here by |x) = 0 x)’ 


N 


riHfQo\x)=ViY. 






-■ p 


N 


-^d^ + Vix,) 


i^A0^)<S)X) 


= Hl^(-l)^P ^ 

P I i=l 

= Vt E(-i)''^ |e ^ [-^d^^PAd^e^ - v’Ad^] 0 x| 


(C2) 


N 


= Pi 1 E b [-di^Ad^e^] 0 X 


At the third equal sign we have used the fact that Vi projects out the wave function belonging to the subspace Qq. 
And at the final equal sign we have used ipAdiO^ = 0, because diO^ generates (5-functions at Xi = Xi±i and the Slater 
determinant ‘fA\xi=xi±i = 0. 

Now let us see how diO^ generates 5-functions. The sector function 9^ can be written into a chain product of step 
functions: 


9^ = 9{x 2 - xi)9{x-i - X 2 ) ■ ■ • 9{xi - Xi-i)9{xi+i - Xi) ■ ■ ■ 9 {xn - xn-i) ■ 


(C3) 
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We therefore have 

diO^ = di [6 {x2 - xi)9{x3 - 3 : 2 ) • • • 9{xi - Xi-i)9{xi+i - Xi)--- 9 {xn - a;Ar-i)] 

= S(xi - Xi-i) [9 {x2 - xi)9{x 3, - X2)--- 9{xi-i - xi-2)9{xi+i - Xi)--- 9 {xn - a;Ar_i)] (C4) 

- S{xi+i - Xi) [9 {x 2 - xi)9{xii - X 2 )--- 9{xi - x^-l)9{xi+2 - Xi+i) - - - 9 {xn - a;Ar-i)] , 
which we rewrite in a simplified notation as 


di9'^ = 6(xi - - S{xi+i - , (C5) 

where 9^-- = 9^I9{xi — Xi-i) is the reduced sector function. 

Let us now consider the summation 5 \—di^Adi9^\. Since di9^ and both generate 5{xi — Xj+i), they 

can be paired up: 


N 


^ [-diipAdi9^] = - 


Z =1 


N 


N-1 


+ Y di<pAS{xi+i - Xi)9Y^^^ 


N-l 


= ^ X! ~ di+ifA)Sixi+i - a;06'[i+i,i] 


(C 6 ) 


N-l 


= Y 9 i^AS{xi+i - Xi) 9 l^^^^^^ 


where, in the last step, we have used ^i+l(pAS{xi-^.l — Xi) = fi+iy [diJ^iLpA^ixi+i — Xi)] = —di^pA5{xi+i — li), where 
£i+i^i is an exchange operator that exchanges the index z + 1 -o- *. Untill now, we have shown that in the identity spatial 
sector a;i < X 2 < • • • < xjq, the operator ViHjQq generates {N — 1) (5-functions of neighboring spatial coordinates. 
The whole expression for PiiL/Qo lx) is then 


ViHfQi,\x)=Y^-^)''P 

P 


'N-l 


Y, diLpAS{xi+i 



(C7) 


Next we act — iLint)~^’Pi on Eq. (IC7I) . Use the fact that when more than two particles are at a same position, 

all the diifiA^s will vanish and so will Eq. (IC7I) . we can deal with the N{N — l)/2 5-functions in Hint and in Eq. (IC7I) 
separately, which means. 


Vi 


E 


(0) 


1 

-Hint 


ViHfQo lx) 


y p 

9 V 


N-l 

Y 

2=1 


5{xi+i - Xj) 
5{xi+i - Xi) 



N-l 


2=1 


diipAdli+i^,-^ 0 X 


(C 8 ) 


In the final step we act (x'| QoHfVi on Eq. (IC 8 I) . (x'| QoHjVi is the hermitian conjugate of a wave function 
having the form of Eq. (ICTl) with a different spin state x! but the same ip a- Look at Eq. ([C7)). since each spatial 
sector has TV — 1 terms, where each term is composed of a 5-function and a reduced sector function, there will be 
totally N\{N — 1) terms appearing in this expression. However, only N\{N— l)/2 terms are of different 5-function and 
reduced sector function. Eor example, consider a sector P (which labels the sector xp-^ < - - - xp^ < xp^_^_-^ - - - < xp„) 
and one of its neighbouring sectors P' = P£i^i+i (which labels the sector xp-^ < ---xp^^^ < xp^--- < xp^^), they 
both possess the term 5{xp.^^ — xp^)9^p_^^ pj. There is also another way to think about this, there are totally 
N{N — l)/2 different 5-functions and for each 5-function there are {N — 1)! different reduced sector functions, so 
totally N{N — l)/2 • {N — 1)! = N\[N — l)/2 different terms. Those different terms are orthogonal to each other, 
because they have different 5-functions and reduced sector functions as well as the fact that when more than two 
particles are at a same position, diipA will vanish. Eor example, one of those N\{N — l)/2 terms belonging to sectors 
P and P' may be 


(-l)^P ^diipA5{xi+i - x{)9Y^^i^ ® x'} + P' [diipA5{x^+i - Xi)9|^p^^^^ ® x'} 

= (-1)^^’ \^^ipA5{x^+l - Xi)9\^p-^ f^ ® x'} - {-P)^P£i,^+l ^diipA5{xi+i - Xi)9Y^^^ ® x'} 

= \diipA5{xi+i - Xi)9\ip-^ i^ 0 [1 - £^,^+l] x'l • 


(C9) 
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Equation (ICSp . similar to {x'\ QoHjVi, also has N\{N — l)/2 orthogonal terms corresponding to different reduced 
sector functions and ‘(5-functions’, as the projection operator Vi plays the role of the (5-functions. So finally, the matrix 
elements of the second-order perturbation effective Hamiltonian Eq. (IB2I) can be evaluated as 

1 


- PiHfQolx) 

^int Hint 






'N-l 


1 ^ t 


^ ^ A^iXi+1 ^i)^[i+l^i] ® [1 £i,i+l] X 




'N-l 


0 [1 - S^^i+l] X 


N-l 


= - - / \\_dxj]rN\ Y \dnpA\^ 5{xi+i - (8){[l-£’*,i+i]x'}'''{[l-£’*,,+i]x} 

d d j=i i=i 

iV-l „ N 

= - -Y HI Y[dxj \d^LpA\^ d{xi+i - a:i)^[Ui.i] ® [1 - X • 

d i=i d 


An effective spin-chain model iJeff is therefore obtained 


Af-l 


where 


HeS =- Y C"* [1 - £i,i+l] , 

9 ^ 


Q = N\ j W_dXj \diipA\^ d{xi+i - . 


(CIO) 


(Cll) 


(C12) 


The above derivation is valid for fermions. In the case of bosons, a general many-body wave function can be written 
as 


T(xi ■ • •XAr,ai • •-CTAr) = 0 x] • 

p 

Following the same procedure as above, we end up with an effective Hamiltonian as 

1 ^-1 

Heff =- Ci [1 + , 

9 ti 


(C13) 


(CM) 


Appendix D: One-body density matrix 

Given a many-body wave function the one-body density matrix is defined as: 

Pa'a{x',x)= Y J dx2---dxN'if*{x',X2---XN,cr',a2---crN)'ifix,X2...XN,cr,a2---<7N)- (Dl) 

<T2---(7N 

For fermionic systems whose wave function takes the form of Eq. ([2]) in the main text, 

d/ = ^x]=^aY ® x] , (D2) 

p p 

where 9^ = P9 is the sector function (generalized step function) for the sector labled by permutation operator P, the 
one-body density matrix can be written as 

Pa'a{x', x)= Y [ dX 2 --- dxN^A^A Y ® {P'x'*){Px) , 

(72---(7]S! P'P 


(D3) 
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where (p'X = ip*^{x' ,X 2 - ■ ■ xn), = ^a{x,X 2 ■ ■ ■ xn), = P'6{x' ,X 2 - ■ ■ xn), 9^ = P9{x, X 2 ■ ■ ■ xn), x'* = 
X*(ct', tT 2 • ■ • ctat): and x = x{^i^ 2 ' ‘ ‘ (^n)- A permutation P can be written as P 2 -N ' where 

is the loop permutation operator defined in Appendix El and P 2 -N is a permutation operator acting on indices 
2,3 ■ ■ ■ N. This means first move particle 1 to position m by a loop permutation, and then permute the remaining 
— 1 particles. Similarly, P' can be written as P' = P^-j^ ■ (1 • • • n). The summation over P' and P can then be 
written into another form: 

Paia{x,x)= ^ f dx2---dxNP’AV’A^ ^ 0 [P2 _n ' {1 ■ ■ ■ m)x*] [P 2 -N ' {I''' n)x] 

P 2 -N 

= ^ f dx 2 ---dx Nip aV’A^ ^ (g) [P 2 -Ar • (1 • • [P 2 -Ar • (1 • • •n)x] 

<T2---<7n rnn P^-n 

= / ^*2 • • • dxNiPAVA9"^'^"''^^9^'^'"'^'> 0 ^ [(1 • • • to ) x '*] [(1 • ■ ■ n)x\ ■ 

mn <72-"<7 n 

(D4) 

The second equal sign follows the fact that if P 2 -N ^ P 2 -N, 9 (x 2 • • • Xm, x' ■ • • xn)9^^-^ {x 2 ■ ■ ■ Xn,x ■ ■ ■ xn) = 0, 

and the third equal sign uses the fact that J dx 2 ■ ■ ■ dxN is invariant under P 2 -n- So the one-body density 

matrix can be separated into a spatial part and a spin part 

Pa'aix , x) — ^ ^ Pm,n{x , x')Sfn,n(i^ j 
m,n 


where the spatial part 


Pm,n{x',x) = {N - 1)1 J dX2---dXN<PA‘PA9'^^"'’^^9^^''^'>, (D 6 ) 

is simply the one-body density matrix of a system of spinless fermions for the spatial sector X 2 < x^ - ■ ■ Xm < x' ■ ■ ■ Xn < 
X - ■ ■ Xn (to < n, for example). And the spin part 

Sm,n{(T',cr) = [(I---to)x*] [(l---n)x] = (x|cL(CT')(”r---n)c„(cr)|x) , (D7) 

(T2---(TN 

where {m - ■ -n) is the loop permutation operator, and cl^(u) can be regarded as fermion (or hard core boson) creation 
operators, which is just a formal symbol to select out the spin states. For bosons, simply change the sectored spinless 
fermionic one-body density matrix to bosonic one by Pm,n = 


Appendix E: Green’s function results for (Nf, Ni) = (1, 1) 

The Hamiltonian of two particles in a one dimensional harmonic trap with a spin dependent magnetic gradient is 

(El) 


1 52 ^ g 2 ^ 2^ 

^ ^ ~ 2 ^ 2 ^^ 2^2+9^X1 -X 2 ) - Gxial - GX 2 CT 2 ■ 


In the absence of the magnetic gradient (i.e., G = 0), there exists an exact solution to the problem [^. Here we 
generalize this solution in the presence of the magnetic gradient. To this end, we make a transformation of operators 
by making spatial and spin coordinates operators into Jacobi coordinates: 


^ a;i - 2:2 


X 1 +X 2 „ CTl - c- ^"1 + 

'? — -:=—, *^2 = 


V2 


V2 ’ 


V2 ■ 


(E2) 


The transformation rules of other operators such as djdx can be obtained from them. The Hamiltonian can be 
separated into the center-of-mass motion part and the relative motion part: 




(E3) 
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For center-of-mass motion, it is a simple harmonic oscillator with center shifted by GS' 2 . For relative motion, it is 
a simple harmonic oscillator with center shifted by GSi plus a (5-function potential at the origin. We can first let 
particle 1 to be spin up and particle 2 to be spin down, then anti-symmetrize the wave function in the end. In this 
case, 5*2 = 0 and S'! = 1/2 are fixed. The eigen wave functions for the center-of-mass motion are still simple harmonic 
oscillator eigen-functions. What matters is the relative motion part. After a coordinate shift Xi — GSi —?► Ai, The 
relative motion Hamiltonian can be written as 

^frei = -i^ + iA2 + ^(5(Ai + G5i), (E4) 

which includes a simple harmonic oscillator part and a (5-function source term. For this relative Hamiltonian, use the 
one-body Green’s function 


G{E-X,,X[) = Y, 

i=0 


1 

E-E, 


ct>,{X,)<j>*{X[), 


(E5) 


where Ei = i -\-1/2 and (pi are the single particle harmonic oscillator eigen-energies and eigen-wave functions, respec¬ 
tively. The corresponding Lippmann-Schwinger equation for the relative wavefunction is given by 


(^(Ai)= J dX,G{E-,X,,X[)^5{Xi+GSiMX[) 

= ^G{E-X,,-GSiM-GSi). 


(E6) 


We just got the expression for the relative wave function, where (p{—GSi) is a constant can be determined by 
normalization of (p{Xi). And the relative energy must satisfy 


G(A;-GS'i,-GS'i) = — . (E7) 

9 

Note that, when G = 0, the Green’s function method fails at E = Ei, and for E ^ Ei the left hand side of Eq. (lE7l) has 
an analytical form . Actually the solution of fully symmetric spin wave function (necessarily assoticated with fully 
anti-symmetric spatial wave function) which has E = Ei should be complemented to the Green’s function solution. 
However, for G ^ 0, there is no such pathological behavior for the Green’s function method. Also to be noted is that 
for one E, there could be only one 1/g for which Eq. (IE7I) is satisfied. This means, for relative motion, there could be 
only one bound state. However, this is no longer true for three particles, because the Lippmann-Schwinger equation 
for three particles is an integral equation and there can exist infinitely many bound states for three particles’ relative 
motion. 

Finally, substitute back Xi —>■ Xi — GSi, after anti-symmetrization, the total wavefunction for two fermions is 
given by 


V2 


^cmiX2) - gV2) in) - ^i-X, - gV2) |;t) 


The center-of-mass separation between the two spins, A = {x\a\ -|- /2, can be calculated as 


(E8) 


A = ^ydAiAiinAOl'+G, (E9) 

where <f{Xi) is decided by Eq. (IE6|) . which dependents on G and E, where E is dependent on G and 1/g by Eq. (IE7I) . 
The first term in Eq. (IE9I) is from interplay between the interaction and the magnetic gradient, while the second term 
in Eq. (|E9ll is due to the harmonic trap shift induced by th e magnetic gradient. _ 
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